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c/j . Abstract 

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma 
waves and their associated instabilities (e.g., Buneman, Weibel and other two-stream 
instabilities) created in collisionless shocks are responsible for particle (electron, 
positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle 
(REMP) code, we have investigated particle acceleration associated with a rela- 
tivistic jet front propagating into an ambient plasma. We find small differences in 
the results for no ambient and modest ambient magnetic fields. Simulations show 
that the Weibel instability created in the collisionless shock front accelerates jet and 
ambient particles both perpendicular and parallel to the jet propagation direction. 
The small scale magnetic field structure generated by the Weibel instability is ap- 
propriate to the generation of "jitter" radiation from deflected electrons (positrons) 
as opposed to synchrotron radiation. The jitter radiation resulting from small scale 
magnetic field structures may be important for understanding the complex time 
structure and spectral evolution observed in gamma-ray bursts or other astrophys- 
ical sources containing relativistic jets and relativistic collisionless shocks. 
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1 Introduction 



Nonthermal radiation observed from astrophysical systems containing rel- 
ativistic jets and shocks, e.g., active galactic nuclei (AGNs), gamma-ray bursts 
(GRBs), and Galactic microquasar systems usually has power-law emission 
spectra. In most of these systems, the emission is thought to be generated 
by accelerated electrons through the synchrotron and/or inverse Compton 
mechanisms. Radiation from these systems is observed in the radio through 
the gamma-ray region. Radiation in optical and higher frequencies typically 
requires particle re-acceleration in order to counter radiative losses. 

Particle-in-cell (PIC) simulations can shed light on the physical mecha- 
nism of particle acceleration that occurs in the complicated dynamics within 
relativistic shocks. Recent PIC simulations using injected relativistic electron- 
ion jets show that acceleration occurs within the downstream jet, rather than 
by the scattering of particles back and forth across the shock as in Fermi ac- 
celeration (Frederiksen et al. 2003, 2004; Hededal et al. 2004; Nishikawa et 
al. 2003, 2004, 2005). Silva et al. (2003) have presented simulations of the 
collision of two inter-penetrating electron-positron plasma shells as a model 
of an astrophysical collisionless shock. In the electron-positron simulations 
performed with counter-streaming jets (Silva et al. 2003), shock dynamics in- 
volving the propagating jet head (where Fermi acceleration may take place) 
was not investigated. In general, these independent simulations have confirmed 
that relativistic jets excite the Weibel instability (Weibel 1959). The Weibel in- 
stability generates current filaments and associated magnetic fields (Medvedev 
and Loeb 1999), and accelerates electrons (Silva et al. 2003; Frederiksen et al. 
2003, 2004; Nishikawa et al. 2003, Hededal et al. 2004). 

In this paper we present new simulation results of particle acceleration 
and magnetic field generation for relativistic electron-positron and electron-ion 
shocks using 3-D relativistic electromagnetic particle-in-cell (REMP) simula- 
tions. In our new simulations, electron-positron and electron-ion relativistic 
jets are injected into electron-positron and electron-ion plasmas in order to 
study the dynamics of a relativistic collisionless shock both with and without 
an initial ambient magnetic field. 

2 Simulation Setup and results 

Four simulations were performed using an 85 x 85 x 320 grid with a total 
of 180 million particles (27 particles/cell/species for the ambient plasma) and 
an electron skin depth, A cc = c/u pc = 9.6A, where co pe = (Aire 2 n c /m c ) 1 ^ 2 is 
the electron plasma frequency and A is the grid size (Nishikawa et al. 2004). 
In two other simulations an electron-positron jet is injected into a magnetized 
and unmagnetized electron-positron ambient plasma and in two simulations 
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an electron-ion jet is injected into a magnetized and unmagnetized electron-ion 
ambient plasma. The choice of parameters and simulations allows comparison 
with previous simulations (Silva et al. 2003; Frederiksen et al. 2003, 2004; 
Hededal et al. 2004; Nishikawa et al. 2003, 2004, 2005). 

The electron number density of the jet is 0.741nb, where is the density 
of ambient (background) electrons. The average jet velocity v- s = 0.9798c, and 
the Lorentz factor is 5 (2.5 MeV). The jets are cold (i>? th = vf th = 0.01c and 
fj th = 0.0022c) in the rest frame of the ambient palsma. Electron-positron 
plasmas have mass ratio m p /m c = m c +/m e ~ = 1 and electron- ion plasmas 
have m\jm c = 20. The electron and ion thermal velocities in the ambient 
plasmas are t>Jr h = 0.1c and v[ h = 0.022c, respectively, where c is the speed 
of light. The time step At = 0.013/u;pe, the ratio u pc /Q c = 11.5, and the 
Alfven speed (for electrons) va c = (^c/^ pc )c = 8.66 x 10~ 2 c. With the speed 
of an Alfven wave given by v A = [Vl/(l + V^/c 2 )] 1/2 = 6.10 x 10~ 2 c where 
Va = [B 2 /Aix(n c m c + ripmp)] 1 / 2 = 6.12 x 10" 2 c, the Alfven Mach number 
Ma = fj/fA = 16.0. With a magnetosonic speed v ms = (t> 2 h + v\) 1 ^ 2 = 0.132c 
the Magnetosonic Mach number M ms = i>j/i> ms = 7.406. At least approxi- 
mately the appropriate relativistic Mach numbers multiply these values by 
the Lorentz factor. Thus, in an MHD approximation we are dealing with a 
high Mach number shock with 7M >> 1. The gyroradius of ambient elec- 
trons and positrons with v± = v t h = 0.1c is 11. 1A = 1.154A CC (for ambient 
ions: 49. 6A = 5.16A ce ). All the Mach numbers with electron-ion jets are ap- 
proximately increased by ^m-Jm c = = 4.47. 

Current filaments resulting from development of the Weibel instability 
behind the jet front are shown in Figs, la and lb at time t = 28.8 /oo vc for 
unmagnetized ambient plasmas. In case (a) an electron-positron jet is injected 
into an electron-positron ambient plasma. In case (b) an electron-ion jet is 
injected into an electron-ion ambient plasma. The maximum values of J y are 
(a) 15.63 and (b) 10.7, respevtively. The electron-positron jet shows larger 
amplitudes than the electron-ion jet at the same simulation time and mag- 
netic fields reduce the maximum values, confirming previous simulation results 
(Nishikawa et al. 2004). The effect of weak ambient magnetic fields affects the 
growth rates of Weibel instability slightly as shown in Hededal and Nishikawa 
(2004). 

The heating and acceleration of jet electrons in directions parallel and 
perpendicular to the flow is shown in Figure 1 for the electron-positron case 
(lc and le) and for the electron-ion case (Id and If). The jet electrons are 
split into two parts: the injected (blue: rear half (Z < 160A) and shocked 
(red: front half (Z > 160A). (The jet electrons are divided at Z ~ 160A).) 
Since in the case of electron-ion jet the Weibel instability grows slightly the 
blue curves in Figs. Id and If are considered as the distributions of injected jet 
electrons. In both parallel and perpendicular distributions, jet electrons are 
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more accelerated in the electron-positron case than in the electron-ion case. 
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Fig. 1. Panels (a, c, e) and (b, d, f) refer to the electron-positron and electron-ion 
cases, respectively (unmagnetized cases). 2D images show the current density (J y ) 
at t = 28.8/u;pe. Colors indicate the y-component of the current density, J y [peak: 
(a) 15.6, (b) 10.7], and the arrows indicate J z and J x . The injected (blue line) and 
shocked (red line) electron distributions are shown as a function of jvn ((c) and (d)) 
and 7fj_ (e and f) where 7 = (1 — (ujj + v'jj/c 2 )^ 1 ^ 2 . 

The electrons are deflected by the transverse magnetic fields (B x , B y ) via 
the Lorentz force: — e(v x B), generated by current filaments (J z ), which 
in turn enhance the transverse magnetic fields (Weibel 1959; Medvedev and 
Loeb 1999). The complicated filamented structures resulting from the electron 
Weibel instability have diameters on the order of the electron skin depth (A ce = 
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9. 6 A). This is in good agreement with the prediction of A m 2 1//4 c7 t 1 / 2 /cj pe rs 
1.188A CC = 11.4A (Medvedev and Loeb 1999). Here, 7 th ~ 1 is a thermal 
Lorentz factor. The filaments are elongated along the direction of the electron- 
ion jets (b) (the z-direction, horizontal in Figure 1). However, in the electron- 
positron jets the current filaments have coalesced in the transverse direction 
in the nonlinear stage. The transverse current (J x ) in the electron-positron 
jets (a) shows significantly more transverse variation than in the electron-ion 
jets (b). 

The acceleration of electrons has been reported in previous work (Silva 
et al. 2003; Frederiksen et al. 2003, 2004; Nishikawa et al. 2003, 2004, 2005; 
Hededal et al. 2004). We see that some of the kinetic energy (parallel velocity 
v\\ « fj) of the jet electrons is transferred to the perpendicular velocity via 
the electric and magnetic fields generated by the Weibel instability as shown 
in Fig. 1. The strongest transverse and parallel acceleration of jet electrons 
accompanies the strongest deceleration of electron flow and occurs between 
z/A — 210 — 240. The transverse acceleration in the electron-positron jets is 
over four times that in the electron-ion simulations. The strongest acceleration 
takes place around the maximum amplitude of perturbations due to the Weibel 
instability at z/A ~ 220 as seen in Figs, la and lb. 

3 Summary and Discussion 

We have performed self-consistent, three-dimensional relativistic particle 
simulations of relativistic electron-positron and electron-ion jets propagating 
into magnetized and unmagnetized electron-positron and electron-ion ambient 
plasmas. The main acceleration of electrons takes place in the region behind 
the shock front. Processes in the relativistic collisionless shock are dominated 
by structures produced by the Weibel instability. This instability is excited in 
the downstream region behind the jet head, where electron density perturba- 
tions lead to the formation of current filaments. The nonuniform electric field 
and magnetic field structures associated with these current filaments deceler- 
ate the jet electrons and positrons, while accelerating the ambient electrons 
and positrons, and accelerating (heating) the jet and ambient electrons and 
positrons in the transverse direction. 

Other simulations with different skin depths and plasma frequencies show 
that the growth rate and spatial structure of current filaments generated by 
the Weibel instability scale with the plasma frequency and the skin depth 
(Nishikawa et al. 2004). An additional simulation in which an electron-ion jet 
is injected into a ambient plasma with perpendicular magnetic field shows 
magnetic reconnection due to the generation of an antiparallel magnetic field 
generated by bending of jet electron trajectories, consequently jet electrons 
are subject to strong non-thermal acceleration (Hededal and Nishikawa 2004). 
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These simulation studies have provided new insights for particle acceler- 
ation and magnetic field generation. Further research is required to develop 
radiation models based on these microscopic processes. 
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